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ABSTRACT 

We present analytical solutions for winds from galaxies with NFW dark matter halo. Wc consider 
winds driven by energy and mass injection from multiple supernovae (SNe), as well as momentum 
- injection due to radiation from a central black hole. We find that the wind dynamics depends on three 

\ velocity scales: (a) ^ (i5/2Af)^/^ describes the effect of starburst activity (and depends mostly on 

■ the initial mass function (IMF) of stars), with E,M as energy and mass injection rate in a central 
^j, region of radius R; (b) v, ^ {GM,/2Ry^^ for the effect of a central black hole of mass M, on gas 

jy^ • at distance R and (c) Vs ~ {GMh/^Crs)^^"^ which is closely related to the circular velocity ivc) of 

NFW halo, with as the halo scale radius and C is a function of halo concentration parameter. Our 
generalized formalism, in which we treat both energy and momentum injection from starbursts and 
radiation from central active galactic nucleus (AGN), allows us to estimate the wind terminal speed 
I— i! to be 2{v1 + 1.5(r — l)v1 — vlY^"^, where F is the ratio of force due to radiation pressure to gravity of 

■ the central black hole. Our dynamical model also predicts the following: (a) winds from quiescent star 
r ^ \ forming galaxies cannot escape from 10^^'^ < Mh < 10^^'^ M© galaxies, (b) circumgalactic gas at large 

, distances from galaxies should be present for galaxies in this mass range, (c) for an escaping wind, the 

■ wind speed in low to intermediate mass galaxies is ^ 400-1000 km/s, consistent with observed X-ray 
temperatures; (d) winds from massive galaxies with AGN at Eddington limit have speeds > 1000 

Q . km/s. We also find that the ratio [2w^ ~ (1 ~ r)w^]/z;^ dictates the amount of gas lost through winds. 

• Used in conjunction with an appropriate relation between M, and Mh, and an appropriate opacity of 

I dust grains in infrared (K band) , this ratio has the attractive property of being minimum at a certain 

Ci • halo mass scale {M^ ~ xg^^"^^-^ M©) that signifies the cross-over of AGN domination in outflow 

properties from starburst activity at lower masses. We find that stellar mass for massive galaxies 

. scales as oc M°-^^, and for low mass galaxies, cx M^^^. 

^ ' Subject headings: galaxies: starburst — galaxies: active — galaxies: evolution 

^: 

^ ! 1. INTRODUCTION 

] In the standard scenario of galaxy formation the baryonic matter falls inside the potential wells created by dark 

si^ ■ matter halos. This in-falling material cools and forms stars. This picture is met with problems as the galaxy stellar 

■ mass function does not follow the halo mass function. Both low and high mass halos have significantly lower than 
' predicted value of stellar masses (Somerville et al. 2008; Moster et al. 2010; Behroozi et al. 2010). To reconcile with 

these problems it has been proposed that the SNe and starbursts provide thermal energy injection and cause large 
amount of mass to flow out of the galaxy as galactic superwinds (Dekel & Silk 1986; Larson 1974; Oppenheimer & Dave 
2006). The star formation is suppressed as the galaxies lose a signiflcant portion of their baryons due to this negative 
feedback. Although this picture can provide explanation for the low ratios of Mi,/Mh for the low mass halos, but for 
^ ■ the high mass halos, gravity becomes strong and the SNe are not sufficient to drive the gas out. In order to resolve the 
discrepancy at high mass end, it has been argued that the AGN outfiows may sweep away the baryons and suppress 
the star formation in high mass galaxies (Silk & Rees 1998; Wyithe & Loeb 2003; Di Matteo et al. 2005; Springe! et al. 
2005; Croton et al. 2006; Bower et al. 2006). These two feedback processes, when considered together, are believed to 
explain the shape of galaxy stellar mass function at both low and high mass end (Binney 2004; Cattaneo et al. 2006; 
Puchwcin & Springel 2012). 

Apart from their cosmological importance as a feedback process, galactic winds have been a topic of research as 
a gas dynamical problem in galactic physics. Speculations on the possibility that galaxies can harbour large scale 
winds followed the models of solar wind developed by Parker (1965). Burke (1968) proposed a model of trans-sonic 
winds from the galaxy with heat and mass addition from SNe. It was further proposed that galactic winds may cause 
ellipticals to lose all of their gas (Johnson & Axford 1971; Mathews & Baker 1971). The review by Holzer & Axford 
(1970) gives an elaborate account of the theoretical aspects of solar and galactic winds. Chevalier & Clegg (1985) 
showed that energy injection at the center can drive a fast super-wind from the dwarf star-burst galaxy M82. In this 
work the gravity of the galaxy was not considered and the obtained solutions were trans-sonic with a heat injection 
up to the sonic point. Wang (1995) modeled the wind from a power law gravitational potential and showed that wind 
may escape the galaxy or settle in a galactic corona depending on the mass of the galaxy and the effect of radiative 
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cooling. Silich et al. (2011) studied the effeets of cooling on winds from individual star clusters with exponential stellar 
density distribution. 

Observations show that the winds do not consist of a homogeneous medium. The hot gas emitting X-rays and 
the cold/warm gas visible in emmision and absorption lines (Strickland et al. 2004; Bouchc et al. 2012; Kornei et al. 
2012), coexist in the galactic winds (Heckman et al. 2000; Veilleux et al. 2005). Hot phase appears as a smooth flow 
of tenuous gas while the colder phase is clumpy in nature. It is usually believed that cold neutral clouds form as a 
result of thermal instabilities in the hot flow and they are entrained with the parent flow because of its ram pressure 
(Heckman et al. 1993). Hydrodynamical simulations with radiative cooling also supported this scenario (Suchkov et al. 
1994; Strickland & Stevens 2000; Cooper et al. 2008). However, recent observations show that the velocity of neutral 
clouds does not correlate with the velocity of hot flow and rather it correlates with the circular speed of the host 
galaxy (Martin 2005; Rupke et al. 2005b). If one considers the momentum driven winds where the radiation from 
the galaxy acts on the dust grains, then these observations can be explained (Murray et al. 2005; Sharma et al. 2011; 
Chattopadhyay et al. 2012; Nath & Silk 2009). However the radiation may not be sufficient in fainter low mass galaxies 
and ram pressure is still required to explain the cold winds in these galaxies (Sharma & Nath 2012; Murray et al. 2011; 
Hopkins et al. 2012). On the other hand in high mass ULIRGs, radiation from AGN may be an alternate mechanism 
for driving outflows apart from SNe and stellar radiation. 

In spite of the general consensus about the AGN driving in quasar outflows, observationally it has been hard to 
establish that outflows in galaxies are also powered by AGN (Rupke et al. 2005a; Westmoquette et al. 2012). However, 
recent observations do show compelling evidence for the AGN driving in galactic outflows (Sturm et al. 2011; Rupke 
& Veilleux 2011; Alexander et al. 2010; Morganti et al. 2007; Dunn et al. 2010; Feruglio et al. 2010; Fu & Stockton 
2009; Villar-Martfn et al. 2011). From the theoretical point of view, considerable amount of work has been carried out 
in modeling radiation driven outflows in the immediate vicinity of the AGN (Murray et al. 1995; Proga et al. 2000; 
Kurosawa & Proga 2009; Risaliti & Elvis 2010). However, the effects of radiation from accreting black holes has not 
been discussed for driving galaxy scale outflows. Murray et al. (2005) proposed the existence of a critical luminosity 
for the AGN (or the galaxy) necessary for the blow-out of all the available gas. The value of this critical luminosity 
depends on the dust scattering opacity at UV. Everett & Murray (2007) studied Parker wind from AGN occurring 
at scales of ~ 100 pc. Debuhr et al. (2012) carried out simulations showing that the initial momentum injection 
and the fast outflow in the vicinity of hole may shock the surrounding ISM and can result in a galaxy scale outflow. 
Whether the outflows due to the AGN are energy conserving or momentum conserving has been a topic of debate 
as well. Silk & Nusser (2010) argued that the energy driven outflows are not possible in the galactic bulges. King 
et al. (2011) proposed that AGN outflows are momentum driven at small scales and energy driven on larger scales. 
Faucher-Giguere & Quataert (2012) showed that the cooling in the region, shocked by AGN radiation pressure, may 
not be effective and the outflows can be energy conserving. McQuillin & McLaughlin (2012) studied the large scale 
motion of momentum-conserving supershells from a dark matter halo. Novak et al. (2012) carried out a radiation 
transfer calculation assessing the efficiency of various components of AGN spectrum in driving outflows. This work 
showed that most of the UV flux is quickly absorbed and re-radiated in IR. 

It is evident from the studies on starburst and AGN driven outflows that these processes play an important role in 
the formation and evolution of galaxies. In galaxy formation models, these two processes are generally invoked using 
simple recipes through feedback factors. However, in theoretical models of winds they have been treated separately 
for the low mass and high mass galaxies. There is a lack of models that develop a complete hydrodynamic theory of 
winds which can envisage both these feedback processes. In the present paper we address this problem analytically 
and bring in both these wind driving agents together and using the analytical results from our calculation we also 
attempt an explanation for the galaxy stellar to halo mass relation. 

We start with the derivation of a general wind equation in §2, which accounts for any possible mass, energy and 
momentum injection. We then present a brief derivation of SNe driven wind model of CC85 and extend this model to 
the case including a dark matter halo in §3. Afterwards in §5 we introduce the momentum injection from the AGN 
and derive a general analytic solution for the galactic wind, which has inputs from SNe injection, NEW gravity and 
the central black hole. This solution leads us to important results like the terminal velocity of winds, condition for 
escaping winds and the dependence of wind properties on the halo mass and the black hole mass. In §5.2 we show the 
velocity, density and temperature of the outflow as a function of the distance from the center. Interestingly a class of 
our solutions can explain the gas reservoirs in the the galactic halos seen in observations (Gupta et al. 2012; Tumlinson 
et al. 2011) and simulations (van de Voort & Schaye 2012; Stinson et al. 2012). In §5.3 we study the cosmological 
implications of our results. We derive the scaling relation between stellar (M^,) and halo mass (Mh) using simple 
recipes and inputs from our models. We compare our analytically derived stellar to halo mass ratio (SHMR) with the 
observations and results of abundance matching. We discuss our results in §6. 

2. BASIC EQUATIONS 

We consider steady spherically symmetric winds driven by mass and energy/momentum injection from processes 
that are confined in a central region of radius R. Consider a heating rate E and mass injection rate M in this region, 
which has a size of a few hundred parsecs (see below) then the basic fluid equations can be written as. 
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Here 1^ is the volume of the central region in which the energy injection and the mass injection is occurring. $(r) 
represents the gravitational potential and f{r) is the momentum injection force per unit mass. These equations do 
not yield a critical point if the heating and mass injection is zero. However, for a finite energy and mass injection 
there is a critical point which is situated right at the boundary of central injection region where the injection turns off 
(Chevalier & Clegg 1985, hereafter CC85). 

In this work we do not consider the radiative cooling, as it is generally believed that radiation loss over the entire 
wind is less dominant than adiabatic losses (Grimes et al. 2009), hence the cooling does not affect the dynamics of the 
wind. However it may still be important for the thermodynamics of the flow (see appendix D) and may result in the 
precipitation of the wind which can not be dealt with the steady stable flow solutions. The dynamics and survival of 
clouds formed by thermal instability in the winds is also an important issue and it has been studied elsewhere in the 
category of cold winds (Murray et al. 2011; Sharma & Nath 2012; Marcolini et al. 2005; Cooper et al. 2009). Therefore, 
we distinguish our model from that of the clumpy winds and in the present work we study analytically the large scale 
dynamics of homogeneous steady wind from a NFW dark matter halo and its cosmological implications. 

By introducing the mach number A4 — u/cs, the above partial differential equations can be transformed to the 
following ordinary differential equation (see Appendix A for a complete derivation). 

M'-l dM^ _2 m mjl+jM'-) ( e/M \ (t + 1) (/(O ' ^) 

M^M^{-f-l) + 2) dr r ^ ^ ' pv 2pv y e{r) J 2(7 - l)e(r) ^' 

22 . . 

where e(r) = + is the specific enthalpy and we will use 7 = 5/3. As mentioned above, the terms E and M (or 

m) represent the energy and mass injection in a central region of size r = R, beyond which they become zero. We use 
R = 200 pc for the present work. 

2.1. Zero gravity case: Chevalier & Clegg's solution 

In this section we briefly reproduce the CC85 solution, therefore this subsection also serves as a consistency check 
for equation 4. In the CC85 solution gravitational force of the galaxy is not considered as the wind speeds were of the 
order of thousand km s~^which is an order of magnitude larger than the circular speed of starburst galaxy M82. There 
is no external driving force (/(r) = 0). Hence the main driving force is the energy injection from SNe in a central 200 
pc region. The size of the central region also marks the boundary where the energy injection and the subsonic part 
ends. The wind equation then can be solved to analytically for subsonic (A^ < 1) and supersonic [M > 1) part of the 
wind. 

In order to derive the subsonic part of the solution, we can use m = ipv/r which results from the integration of 
continuity equation. Therefore equation 4 becomes 

_i dM"^ _2 3(1 + 7X2) ?,{l + -iM^) ( E/M \ 

1 (5) 



M'^{M^{-1 -I) + 2) dr r r 2r \ e(r) 

Direct integration of energy equation 3 by retaining the energy and mass injection and neglecting the gravity and 



external driving results in, e(r) = E / M = 2v1 , where we have defined a velocity parameter, w^, = y E /2Ad . It can be 
shown that at the critical point, Vcnt ~ '^(R) = Cs{R) = i'^. Substituting this e(r) in equation 5, we get, 

- 1 dM^ _ -1 - 
7W2[(7- 1)7W2^_2] dr r ^' 

This can be integrated to get the following solution, 

^<(^)=i^7T^j [ 7 + 1 ) =Jl-^ for r<R (7) 

Following similar steps, we additionally drop the injection terms to evaluate supersonic part of the solution. Ana- 
lytically it means that we should set $ = / = M = £' = Oin the wind equation 4. Therefore, we are left with the 
following differential equation, 

dM^ 2 

(8) 
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This can be integrated to get the fohowing solution for the supersonic part of the wind, 

These expressions were arrived at by CC85. We shah use these definitions of 5^{Ai) and (5>(A^) in the rest of the 
paper. The terminal velocity of the wind in this solution can be obtained easily by the energy equation which is 
e(r) = + Cs/(7 — 1) = 2w^. For r — > oo, the sound speed can be neglected and one gets Voo = = (2£^/M)^/^ . 
In the supersonic part of the wind for M > 1 the relation in equation 9 can be approximated hy oc r^. Therefore 
when the velocity attains its terminal value then cx T cx r^'^^^ and p cx r~^. 

3. SN DRIVEN WINDS FROM NFW HALO 

In this section we study the effects of NFW dark matter halo on the wind velocity in detail. One can therefore 
consider the calculation in this section as an extension of the supersonic part of CC85 solution. We proceed first by 
fixing the dark matter halo parameters in the next subsection. 

3.1. Dark matter halo properties 

Here we discuss the properties of the NFW dark matter halo which will be used throughout the rest of the paper. 
We consider a Navarro Frank and White(NFW) dark matter halo with a density profile p = psrf /r(r + r^)^ (Navarro 
et al. 1996). For a dark matter halo of total mass Mh, which collapses at a redshift z, the virial radius is given by 



r. 



vir 



[3Mh/i4.7:Ac{z)n^^Parit)]'^^^{l + z)-'^ Where is the critical over density which can be written as Ac( 



IStt^ — 39(i^ + 88dz, with dz = fl^ — 1 and fl^^ = ^ — ( i ^. ^ yi'+nt+o ( i + z ) ^ with parameters having their usual meaning 
(Bryan & Norman 1998; Bullock et al. 2001). We set fl^n = 0.258, fl^ = 0.742, h = 0.72, according to the results 
from WMAP5 (Komatsu et al. 2009). Using these values, the virial radius can be written as 
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The corresponding circular speed is, Vc = \/ GMh/r^i-^. So it can be written as, 

1/6 

143 ( ::" rl (l + z)i/2 ]^ms-i (11) 
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The gravitational potential due to the NFW dark matter halo can be written as, 

GMh ln(l + r/r,) ^ Ml + r/r,) 

^ "ln(l + c)-c/(l + c) r = "^"^ r/r. ^^^^ 

In this v'l = GMhjlCr ^ and C — ln(l + c) — c/(l + c). Also — r^-^jc is the NFW scale radius with c as the halo 
concentration parameter which depends on both halo mass and the redshift of virialization. We use the fitting formula 
for c given by Muiioz-Cuartas et al. (2011). 

3.2. Effect of dark matter halo on winds 

We study the effect of gravity in the supersonic part of the solution. In the subsonic part the NFW gravity does 
not make a difference as shown in Appendix B. Therefore the subsonic part of the previous section is still valid and 
wc can assume that the energy injection propels the gas to a speed of Ucrit = = {E/2MY/'^ at the critical point. 

To study the supersonic solution beyond the critical point {r > R) , we set E and m equal to zero, and use the NFW 
potential in equation 4, which results in the following equation for the supersonic part, 

M'^-l dM' _2 2 rf$NFw(r) ^^g^ 



M^[i-f-l)M^ + 2] dr r e(r) dr 
where e(r) = /2 + Cg/(7 — 1). From the direct integration of energy equation we obtain 

^,2 ^2 



*NFw(?') 



= 

R 



2 7-1 

^e(r) = 2w2 + «'NFw(i?)-*NFw(f) (14) 

where we have used the value of velocity at the critical point as Ucrit = v{R) ~ Cs{R) = and Vi, ~ {E/2M) 
Substitution of this e(r) in equation 13 followed by integration results in, 

\n\5>{M)\ = 21n|r| + 2 In |2u2 + $nfw(-R) - $nfw (0 1 + consi. (15) 
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where i5>(A^) is defined in equation 9. In Appendix C we have plotted the contours of M versus r for an arbitrary 
value of constant. For the analysis here we pick the solution which has a critical point at radius R. To pick that 
solution we set, Sy{M) = 1 and r = i? in the above equation which fixes the constant and we get, 
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Using the definition of $nfWj in equation 12 we get. 
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where we have defined a velocity scale Vs = GMh /2Crs , which depends on the halo mass Mji and the concentration 
parameter c. Comparing with eqn 9, we can clearly see that the solution picks up an additional term due to the ef- 
fect of NFW halo. The solution has three different regimes depending on the value of the ratio v1/v1, as described below 



• case 1 : Vi, » Vg 

This is the case when the initial injection upto the critical point is very strong and hence the wind velocity at 
the critical point is large. Then we can neglect the ratio w^/uj, and we get 

Sy{M) = {r/R)^ (18) 

The NFW halo has negligible effect in this particular case and we have recovered the supersonic part of CC85 
solution. 



• case 2 : v^, = Vs 
In this case we get 

6>{M) = {^ ln{l + r/r,)y (19) 

which shows that the gravity of dark matter halo affects the flow significantly. In the large M limit we find that 
the mach number and the velocity increases slowly with r, as compared to the case without gravity, because of 
the logarithmic dependence in the above equation. 



• case 3 : v^, < Vs 

In this interesting case, we find that the wind speed decreases and finally becomes zero at some distance. This 
distance is decided by the ratio v1/v1, and can be determined by requiring that R.H.S = at some r = Vf. We 
get, 

M1±I£M^1_^ (20) 

If the ratio v1/v1 ^ then Vf ^ Q from the above equation, and hence the velocity is zero at the base of the 
wind, which contradicts the assumption in deriving this solution that fiow becomes supersonic. It means the 
ratio needs to have a finite value so that we get avf > R and in this case velocity starts from being equal to the 
sound speed at r = R and then stops at a distance r/. 

Another interesting result which can be deduced from the above equation is the condition for escaping wind. If 
Vn = Vs then rj = oo, which means that winds will escape the galaxy. Hence any heating occurring at the base 
need to somehow propel the gas at least to a speed of ~ Vg so that the gas can escape the galaxy. 

In figure 1 we plot the supersonic wind velocity for a reference halo mass Mh — 10^^ Mq with other NFW parameters 
calculated using definitions in §3.1 at z = 0. The value of Vs = 236 km s~^for this halo. The solid line is for v'^ = IQv^ 
and in this case the injection is large so that the gravity of the halo does not affect the flow. That is why the solid 
line leaves the plot window and the wind velocity at largest distance is ~ 1500 km s~^which is double the value at 
r = i? in agreement with the result Voo = 2vn of previous section on the solution of CC85. The dashed line represents 
the case 2 when Vi, = Vs and the NFW gravity is important. Consequently the velocities are reduced and at the final 
distance we get a wind velocity ^ 200 km s^^. The dotted line is for = O.^v^ and it represents the case 3 discussed 
above. This line exhibits the flow which ends inside the galactic halo at a final radius, r/ ^ 50 kpc. We find that for 
the low values, which implies a low value of velocity at the critical point, the gas stops at some point within the 
halo. We note that similar solutions were obtained by Wang (1995) for power-law potential, and for specific values 
of wind speed and mass loss rate, but the calculations were done numerically. The analytical model presented here 
generalizes these solutions to NFW halo and to general values of wind speed. 
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Fig. 1. — The wind velocity in the supersonic regime for the solution with NFW halo. The reference mass used is Mh = lO^^M©. The 
solid, dashed and dotted line are for v'^/v'l = 10, 1, 0.3 respectively. The x-axis starts at r = R = 200 pc and ends at the virial radius of 
the halo. 



4. INJECTION PARAMETERS AND VALUE OF 

In the last section we discussed the importance of the ratio, v1/v1 = E / {2Mv1). Here we would like to determine the 
possible values of Vi, which we will use for the rest of the analysis. We can define the energy injection from supernovae 
per unit time as Esn = ch/snM-i, x (10^^ erg) where 10^1 erg is the energy output of a single supernova, a is the 
fraction of this energy retained by the gas after radiative energy losses, M^, is the SFR in solar mass per year and fsN 
stands for energy injection per year per solar mass of star formation. For a Kroupa-Chabrier IMF fsN ^ 1-26 x 10~^. 
Therefore we have, 

erg s-i (21) 
The mass injection rate is written as 

M = (3TZfM^ (22) 

Where TZf is the return-fraction. Typically 30% of the mass is returned to the ISM hence 72./ w 0.3. The factor (3 
takes into account the entrainmcnt of mass from ISM which can increase the overall mass injection. For the starburst 
galaxy M82, /3 is in the range 1.0 to 2.8 which gives M ^ 1.4 - 3.6 Mq yr'^ (Strickland & Heckman 2009). Martin 
(1999) found that for a range of galaxies from low to high masses the mass loss in hot wind roughly scales with SFR . 
We therefore use M ^ M-^, which corresponds to /3TZf ^ 1 . 

With these value of Esn and M we estimate the value of which is {Esn /'^MY^'^ — 562-^^ km s~^. Wc consider 
two modes of energy injection from SNe. In the first case, almost 90% of the energy of SNe is lost via radiation and 
only a small fraction goes into heating of the wind. For this 'quiescent mode', we use a = 0.1 which gives sa 180 km 
s^^. In the other case, when the central injection region is dense and the supernovae are overlapping so that radiative 
losses reduce, and due to thermalization, 30% to 100% of the supernova energy goes in heating the wind (Strickland 
& Heckman 2009). This type of situation is evident in galaxies like M82 where the SFR is generally high and the 
injection regions are supposedly quite dense. For such a case a = 0.3 - 1.0 and we get « 308 - 562 km s^^. To 
represent this mode we take = 500 km s~^. We will call this the 'starburst mode'. 

We would like to point out that the values of Vi, chosen here brackets the range of the possible values. We emphasize 
that = Cs ~ 500 km s~^ corresponds to a temperature ^ 0.7 keV, which is consistent with the hot wind temperatures 
observed in a wide range of galaxies (Martin 2005). On the other hand the quiescent mode (w* = 180 km s~^) is suitable 
for the galaxies with low values of SFR, like our own Galaxy. This particular mode also yields interesting results, and 
we will discuss the implications in a later section. 

5. WINDS IN THE PRESENCE OF AGN 

The AGN is important as it gives a strong momentum injection to the gas via its radiation field. A large fraction of 
AGNs show the evidence of outflowing gas, and it is possible that all AGNs drive outflows and they are observed when 
they are viewed edge-on. Theoretically, these outflows have been associated with the co-evolution of black holes and 
the bulge of the host galaxy (Silk & Rees 1998; King 2003, 2005). How the AGN interacts with the ISM of the host 
galaxies and whether it can drive a large scale outflow which can escape the galaxy is an important question. AGN 
can affect the gas in the host galaxy indirectly where it produces fast nuclear winds which shock the ISM into shells. 
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The fate of these shells is then decided by the supply of energy and momentum injection from the inner regions. Apart 
from this indirect way AGN can also interact with the dust-rich ISM directly via its radiation field. This interaction is 
capable of driving large scale outflows (Murray et al. 2005). Here we consider this mechanism and model the outflows 
as being driven by (Eddington limited) continuum radiation from the black hole. 



5.1. Effect of momentum injection from AGN 

Momentum injection can be provided by the AGN in several ways. Firstly it can be provided via the scattering of 
photons by the free-electrons. As the Thompson opacity is generally small {kt ~ 0.34 cm^ gm~^ for a fully ionized 
gas with solar metallicity), this may be effective in the regions close to AGN where the densities are quite large and 
the radiation field is strong. Another way of momentum injection is to produce outflows via line driving mechanism 
as discussed in (Murray et al. 1995). We consider the momentum transfer to dust grains by the AGN radiation. Dust 
opacities are rather high and recent models of momentum driving of outflows due to AGNs and galactic radiation as 
well, consider the scattering of photons by dust grains. The gas in turn is coupled to these charged grains and gets 
dragged. 

Let us derive the force due momentum injection by AGN radiation (/(t)) which then will be substituted in equation 
4. In its general form The momentum injection from a radiation field can be written as, 

/(r) = K- (23) 
c 

where k is the volume averaged opacity, J- stands for the frequency integrated flux of radiation and c is the speed of 



2 



light. For spherical symmetry and for a optically thin atmosphere, the radiation flux can be written as = L/Anr 
Hence the force of radiation becomes KL/Airr^c. This force has an inverse square dependence on r, hence it can be 
represented as a factor (F) times the gravitational force of the black hole. Therefore, if the gravity of the central 
black hole is given by fgrav = —GM,/r'^ then in the presence of an outward radiation force, the effective force is 
written as fg,eff = — (1 — F)GM,/r^ with F = k,L/{AttGM,c). For F = 1, the effective force is zero, and for the 
case of Thompson scattering the corresponding luminosity is called the Eddington luminosity (Le)- The luminosity 
required to exactly cancel the gravitational force may be different depending on the opacity and process responsible 
for momentum injection. For example, in case of scattering of UV light by dust, for which the opacity is roughly 
3500 times the electron-scattering opacity (Drainc 2003), only a luminosity of ^ 0.001 Le is required to counter the 
black hole gravity. However, it has been showed recently that most of the UV photon field from the AGN may get 
attenuated within a short distance because of the large optical depths in AGN environments (Novak et al. 2012). In 
that case the re-radiated Infrared(IR) photons serve as the mainstay of AGN radiation (see also Dorodnitsyn et al. 
(2011). Therefore in this work we consider the momentum injection from IR radiation. IR to dust scattering opacity 
is, Kjji = 13 cm^ gni~^, in K band (Li & Draine 2001; Draine 2003) for a gas to dust ratio of 125. Opacity in IR is 
not as large as it is in UV, however we find that it is large enough to drive strong outfiows in massive galaxies. Hence 
the momentum injection force in our case becomes, 

m -^ = V^ (24) 

where F = kjuL/ {AttGM,c). We consider an Eddington limited AGN [L = Le), where Le = {4:TtGM,c/ kt) is the 
usual Eddington luminosity for electron scattering. This fixes the value of F = kir/ hlt ~ 38 , for our case. 

We need to justify that the atmosphere is optically thin in IR so that we can work with a constant F. In order to 
do this we estimate the optical depth for IR light. Optical depth can be estimated as, r ^ J Kp{r)dr. Using Km — 13 
cm^ gni"^ and p ^ 10 mp cc^^which, as we will see is an upper limit for density in our wind models, the optical 
depth comes out to be 0.01 at 20 pc and 0.13 at 200 pc. Even at the edge of the injection region the optical depth is 
very small, hence we can conclude that the atmosphere is optically thin to IR radiation. Beyond 200 pc, the density 
decreases rapidly, as the wind expands adiabatically in the supersonic regime. As we will see in next section, density 
becomes as low as 0.001 mp cc~^ at 10 kpc, hence the wind material stays optically thin to IR radiation at large 
distances as well. 

Let us consider the supersonic section of the wind in which both SNe and AGN radiation are effective. (See Appendix 
2 for the subsonic part of this wind.) For the supersonic part the energy and the mass injection both are zero, although 
the gravity due to NFW halo and the effective force due to radiation and gravity from the central AGN are present. 
The total potential can be written as, 

$,„,,,(r) = ci>.(r) + ci>NFw(r) = - 2v'^ ln{l + r/rs) ^2^^ 

r r/rs 

With the aid of this total potential and the momentum injection term /(r) TGM,/r^ , the wind equation 4 for 
supersonic part (r > R) can be written as, 

2 2 ^(l-F)GAf. , d^NFw{r)\ 
- ' (26) 



M'^[{'^ -l)M'^ + 2] dr r e(r) V dr 
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Next wc integrate the energy equation directly as below, 

r 

= 



, , (1 - V)GM, , ^ , . 
e(r) - ^ + ^NFw{r) 



R 

e(r) = [2vl^, - 2(1 - T)vl + $NFw(i?)] + 2(1 - T)vl^ - $nfw(0 (27) 

where we have used e(i?) = 2u^j,;j = 2v^ " (1 ~ r)v^, from Appendix B for the subsonic section of this wind. We have 
also used = GM,/2R. The term inside the square bracket in the above equation is a constant, so if we substitute 
eq 27 in 26 it becomes exactly integrable and we get, 



\n\5y{M)\ = 21n|r| +21n 



2vl,, - 2(1 - T)vl + $NFw(i?) + 2(1 - T)vl^ - $NFw(r) 



const. (28) 



Using the condition^, 5y{M) = 1 at r = i?, and substituting the expression for NFW potential we get. 



5>{M) = [r/Rf 



1 _ (1 _ T)4- fl - :^ V — f Mi±i^ _ H^ + r/Ts) 



(29) 



This equation gives the complete solution of the wind from a galaxy driven by SNe and AGN. Let us discuss a few 
asymptotic bahaviours, 

Terminal speed : 

We use the Bernoulli equation 27 to obtain the terminal speed. Taking the limit r — > oo in equation 27 and 
neglecting the sound speed we obtain the following general expression for terminal speed of SNe and AGN driven 
wind from a NFW halo, 

o X 1/2 

,2 I '-*/T^ 1^„,2 „,2 



Voo^2^,t + -{T-l)vt-vij (30) 

Here v1 = GMh/2Crs. In the absence of dark matter halo and the AGN we can neglect Vg and f, which gives 

Voo = 2wJ. However in presence of NFW gravity but no AGN we get the relation Vao = 2\/tJ— wj- For dwarf 
galaxies, the effect of NFW gravity can be neglected, and therefore the wind speed in starburst mode is expected 
to be Voo ~ 1000 km s~^, consistent with observations of winds in starburst galaxy M82. If the black hole is 
massive so that the v, dominates then by neglecting and Vs and putting F = 38, we get v^o ~ 15u,, which for 
a black hole of mass ^ 10^ Mq estimates to Voo ^ 1500 km s^^. 

• Behaviour with r : 

In equation 29 if the SNe injection term (vi,) is dominant then we get (x r^, T (x r~'^^^ and p cx as 
derived in §2.1. On the other hand if the AGN provides the main driving then neglecting Vg and v^, from equation 
29 we get, Ai^ oc r. In that case when the velocity attains its terminal value, we obtain oc T cx r~^/'^ and 
p cx r~^. 

It is clear that, apart from the f^, the wind properties depend on two more terms. One involves v, which is a 
parameter of the black hole and the other being Vs, which depends on the mass of NFW halo. To proceed further we 
need to learn whether the black hole mass is somehow related to dark matter halo. 

Observations show that the black hole mass is related to the velocity dispersion of the central spheroidal bulge 
component of the galaxies, and the relation can be written as 

''>^'''='^''<200^) ^''^ 

This correlation has been extensively studied in the literature, with slightly different values for a and b. We will list 
a few of them. Gebhardt et al. (2000) was one of the earliest to report the correlation with a ~ 8.08 ± 0.07 and 
b = 3.75 ± 0.3. Ferrarese (2002) gives a ~ 8.22 ± 0.08 and b = 4.58 ± 0.52. Values in Tremaine et al. (2002) reads 
a = 8.13 ± 0.06 and b = 4.02 ± 0.32. For our analysis we use a = 8.12 and b = 4.24 from a recent study by Giiltekin 
et al. (2009). 

We also use a relation between the a and circular speed Vc- Ferrarese (2002) reported a correlation between the 
circular speed of the galaxy (wc) and the cr, given roughly as Vc cx u^-^'^^om^ Similar relations were also deduced by 
Baes et al. (2003) and Pizzella et al. (2005). All these relation are very close to the linear relation Vc = ^/2a for 
spherical mass distribution (Binney & Tremaine 2008), found in massive ellipsoidal galaxies where the bulge to total 

mass ratio is unity (see also Volonteri & Stark 2011). Hence for our analysis we use Vc 

The relation between Vc and a breaks down for low mass galaxies (Ferrarese 2002). It is easy to understand this as 
the lower mass galaxies may admit a bulge to total mass ratio less than unity. Hence the galaxies with smaller bulges 

^ By imposing this condition wc have picked the solution which becomes supersonic at r = R = 200 pc. To see the complete solution 
space with a variable constant the reader is referred to the M versus r diagrams in Appendix C. 
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will have the mass of their black holes smaller than the one expected from M^-a-Vc relation. On similar grounds, 
Kormendy & Bender (2011) concluded against the co-evolution of central black hole and dark matter halo. However, 
they also showed that a cosmic conspiracy causes Vc to correlate with sigma for massive galaxies {vc ^ 200km s~^). 

We would like to emphasize here that even if the v^-cr relation breaks down below a Uc ^ 200 km s"-'^, still it does 
not make a difference in our results because the AGN term is effective only for very massive systems. In equation 
26 it is the relative values of v, and Vg, which govern the dynamics. For Uc=200 km s~^, using c = 12 we get 
Vs = {GMh/ICrsY^"^ ~ 400 km s^^which is larger than {^/T — l)v, 100 km s^^. and the black hole term is even 
lower for < 200 km s"^. 

Using the Af.-cr relation provided by Giiltekin et al. (2009), aided with the relation a = Vc/V2, where v^. is given 
by equation 11 we arrive at the following relation between black hole mass and the halo mass (see also Volo uteri & 
Stark (2011)). 



M. ( Mu 



0.74 X 107 MqJ VlO^^ Mq yr 



187r2 



1/2 



(1 + Zvir)^/' (32) 



where Zvir is the redshift at which the halo collapsed and got viriallized. This relation can be used to determine w, 
as a function of halo mass through, w, = {GAI,/2R)^^^ . Substituting v, and = \/ GMhj^r JZ in equation 29 with 
Ts and C determined from the redshift dependent definitions for NFW parameters given in §3.1. This enables us to 
compute the mach numbers as a function of r for the galaxy of any desired halo mass(Af?i) and its halo collapsing 
at any desired redshift. The mach numbers can then be converted to the velocity simply by multiplying it with the 
sound speed which can be obtained using the relation = e(r)/(AlV2 + 1/(7 - 1))- 
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Fig. 2. — Wind velocity at virial radius as a function of halo mass for four different redshifts, z = 0, 2, 4, 6 in the panels from top- left to 
bottom-right respectively. The thin red lines in each panel are for starburst mode (i;* = 500 km s~^) and thick blue ones are for quiescent 
mode (ti* = 180 km s^^). Dashed line represent the outflow speeds without AGN. The green bar on x-axis in top two panels gives the 
range of galaxies in which the gas reservoirs in halos are observed (Tumlinson et al. 2011). 

In figure 2 we plot the wind velocities at virial radius^ (wwind,vir) as a function of halo mass(Mh). The velocities are 
obtained by solving equation (29) with inputs from equation 32 and definitions in section 3.1. We show the results 
corresponding to four different value of ^vir iu four panels, red lines are for the starburst mode = 500 km s~^and 
the blue lines are for the quiescent mode = 180 km s~^. The dashed lines represent the solution without AGN and 
these show that the wind velocities cut of at some halo mass as the gravity becomes strong enough to counter the 
energy injection. For larger value of this cut-off occurs at a larger halo mass. 

^ Wind velocity at virial radius (^'wind vir) need not be equal to the the terminal speed at infinity (uoo) as the later is calculated by using 
the fact that the gravitational force at the infinity is zero while in case of «wind vir there may be a contribution from NFW gravity at virial 
radius. 
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Consider a situation in which for low halo masses the AGN driving term is smaller than the NFW gravity term. 
This implies that the wind velocity at virial radius in low mass halos decreases with increase in halo mass. However 
the black hole mass also increases with the halo mass. Since the slope of M,-Mh relation is greater than unity, hence 
the rate of increase of the black hole or AGN term is larger and at a particular halo mass it overcomes the NFW 
gravity term. Hence for largest galaxies we should see an increase in wind velocity with halo mass, which indeed is 
the case as shown by rising solid lines for massive halos. One can further compare the solid lines with dashed lines 
where the latter represent the case without the AGN and does not show any winds in high mass galaxies, which is a 
confirmation of the fact that in high mass galaxies outflows are driven by AGN. 

The thick solid blue curve is special as it features the wind quenching due to NFW gravity as well as the high 
velocity winds due to AGN. There is a falling part of the curve which exhibits the cut off at some halo mass and then 
there is another part which rises at some larger halo mass. Hence there exist a range of halo masses roughly within 
10^^~^^Mq, which do not have escaping winds. The exact value of this range depends on the redshift of collapse (zvir)- 
For example at z^ii = the rising part of the thick blue curve, which shows the effect of AGN, starts rising beyond 
a halo mass of lO^^M© while for z^ir = 2 it rises roughly at a halo mass ~ IO^^Mq. This is easy to understand as 
the AGN driving depends on the black hole mass which does increase with redshift (see equation 32). Also the falling 
part of the thick blue curve ends at a smaller halo mass for a larger value of Zvir, because the value of increases 
with redshift. We would like to mention here that the recent detection of gas reservoirs in the halos of galaxies by 
Tumlinson et al. (2011) covers roughly the similar range in halo masses shown by the green bar on x-axis in upper two 
panels of figure 2. 

If we focus on the wind velocities in lower to intermediate mass halos Mh < IO^^Mq, we find that AGN never 
dominate in these and if there are winds they have to be driven by starbursts and SNe. We note that the wind 
velocities in low mass galaxies fall in the range 400-1000 km s~^, depending on the efficiency of the energy injection 
process. These velocities agree with the ones inferred from the X-ray temperatures of the superwind regions in dwarf 
straburst and luminous infrared galaxies (Heckman et al. 2000; Martin 1999). However in case of galaxies with halo 
mass, Mh > IO^^'^Mq the wind velocities either exceed 1000 km s~^or they are quenched, depending on whether the 
AGN is present or not. 

We note that in low mass galaxies where the outffows are driven by SNe, the wind velocity, Uwind ^ 1000 km s^^. 
However this limit is exceeded when an AGN is present, since the curves with AGN show wind velocities > 1000 km 
s~^. This reveals the presence of a dividing line of 1000 km s~ ^between SNe and AGN domination in velocity space 
as well. 

In a hot wind with velocity of Wwind the neutral clumps in the wind can be dragged via the ram pressure. Maximum 
velocity these clouds can achieve is the velocity of the hot wind. As mentioned above, Uwind is always within 1000 
km s~^at the low mass end where SNe injection dominates, therefore the cold clouds should also be outffowing with 
velocity lower than 1000 km s~^. On the other hand at the higher mass end, where the AGN dominates, the velocities 
may exceed this limit. Interestingly, the observed outflow speeds of the neutral clouds in diverse galaxies like dwarf 
starbursts, ULIRGs, post-starburst galaxies and Low-ionization BAL quasars also follow this trend (Tremonti et al. 
2007; Trump et al. 2006). 

It has been debated in the literature that the neutral cold/warm winds in ULIRGs are driven by ram pressure of the 
hot wind and/or radiation from stars in the galaxy or by the AGN. If we consider the winds driven by stellar radiation 
then the wind speed is roughly 3 times the circular speed of the galaxy (Murray et al. 2005; Martin 2005; Sharma 
et al. 2011). For a massive ULIRG with a circular speed ~ 300 km s~^the wind velocity predicted by radiation driven 
wind model will be 900 km s~^. On the other hand, if ram pressure is the driving mechanism, then also the velocities 
of the smooth wind and hence of the neutral clouds can not exceed 1000 km s^^unless an AGN is present. Therefore 
we emphasize that the observations of wind velocity in excess of 1000 km s^^indicate the presence of an AGN. 

5.2. Wind properties with distance : Implications for gas observed in galactic halos 

In last section we have established that for a particular mass range the galactic winds may not escape the galaxy. 
Therefore these galaxies are not very important for the intergalactic medium (IGM) enrichment. Interestingly our 
Milky Way with a total mass roughly ~ IO^^Mq, also falls in this mass range. However an important question arises 
for these type of galaxies, as to whether or not these galaxies can retain all the gas in the disk even though they can 
contain the gas inside rvir- 

We show in this section that infact these galaxies have outflows which spill their gas reservoir throughout the halo. 
In figure 3 we show the radial dependence of wind properties for three galaxies which differ in the values of their halo 
mass. We have considered the halos collapsing at a redshift of z^ir = 2 which corresponds to a look-back time of 
roughly ~ 10 Gyr. We choose this collapse redshift, as a fiducial value, in order to model galaxies like Milky Way and 
more massive galaxies, whose halos were already in place by z ^ 2. The thick disk in our Galaxy shows that Milky 
Way underwent its last major merger before ~ 10 Gyr, which corresponds to z ^ 2 (Gilmore et al. 2002). Therefore 
our results can be compared with observations of winds at low redshift {z < 2) universe. 

In figure 3, the thick blue lines refer to the quiescent mode (w* = 180 km s^^) and the thin red lines represent 
the starburst mode (w* = 500 km s~^). The solid lines denote the effect of AGN activity and dashed lines consider 
only the SNe injection. The upper panels plot the velocity as a function of radius, and the middle and bottom panels 
plot the temperature and densities respectively. The densities are estimated using the relation p = M/iirvr'^, where 
M ~ Mi,. For three galaxies with Mh = 10^^, 10^^, lO^^M©, we have used M^ = 1,3, 10 Mqji'^ for quiescent mode 
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Fig. 3. — The wind properties as a function of distance from the base. Legend is same as in figure 2. From left to right, the three columns 
correspond to three different halo masses, representing a dwarf galaxy, a Milky Way size galaxy and a giant galaxy respectively. The x-axis 
in all the plots extend upto the corresponding virial radius. The solid lines represent the winds in presence of an Eddington limited AGN 
and the dashed lines, the case of no AGN. Thin red lines represent the starburst mode (ri* = 500 km s~^) and the thick blue lines are 
for the quiescent mode [Vi, = 180 km s~^). The middle row contains the temperature profiles of winds in three galaxies. The lowermost 
row shows the densities. To calculate densities we have used M^, = 1, 3, 10 M© yr^^ for v^, = 180 km s"\ and = 10, 30, 100 Mq for 
V* = 500 km s~^, for three galaxies respectively. 

and Af,t = 10, 30, 100 MQyr~^ for starburst mode. 

The curves show that for low mass galaxies, all types of winds (with or without AGN, quiescent and starburst mode) 
escape the virial radius. At the other extreme, for massive galaxies, winds with AGN activity can escape, but without 
an AGN, they stop at a distance within the halo (~ 10-200 kpc). The gas temperature and density slightly rises at 
this final halting point due to adiabatic compression. We would like to mention here that the radiative cooling can be 
important for these particular cases as the steady solution ceases to exist beyond a point. Even if the cooling time is 
shorter than the flow time initially, after many flow crossing times the cooling will become effective which may cause 
thermal instability. This can lead to formation of clouds whose fate then will be decided by the physical properties in 
their environment (see also Wang 1995). 

For the intermediate mass galaxy {Mh ^ IO^^Mq), we find that both of the thick blue lines (i.e., quiescent mode of 
star formation with or without AGN) arc contained inside the virial radius. This implies that the wind needs strong 
starburst activity in order to escape the galaxy irrespective of whether or not AGN is present. The dashed thick blue 
line corresponds to a quiescent mode of star formation without AGN and roughly corresponds to our own Galaxy. 
Interestingly, we find that a slow wind is possible, which may extend to a distance of ^ 20 kpc. This can explain the 
recent observations of clouds roughly at 10-20 kpc in our halo (Keeney et al. 2006). 

The solutions which end inside rvir are important in the wake of recent observation of warm-hot gas clouds in the 
halo of our galaxy (Gupta et al. 2012) and for other galaxies as well (Tumlinson et al. 2011). Also recent simulations 
confirm these gas reservoirs around the galaxies of intermediate masses. Recent works find that although in the general 
scenario of galaxy formation the intermediate mass galaxies are efficient in retaining their baryons, these galaxies do 
not retain all of the baryons. It appears that only 20% to 30% of baryons are converted to stars in these galaxies as 
well (Somerville et al. 2008; Moster et al. 2010; Behroozi et al. 2010). Our results provide a natural explanation for 
the missing baryons in these intermediate mass galaxies. 



5.3. Cosmological implications 

Here we derive a relation between the stellar mass and the halo mass using the results for SNe and AGN driven 
outflows in the previous section. In the the scenario of hierarchical structure formation low mass galaxies form at 
earlier times and post formation history is influenced by mergers and the periods of enhanced star formation activity. 
The semi- analytical modeling (SAM) (sec Baugh 2006), which uses simple recipes for feedback and follows the structure 
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formation according to ACDM model, can explain the observed properties of galaxies (Somerville et al. 2008). In recent 
years their have been a growing amount of observational evidence that the massive black holes were already in place 
at high rcdshifts. Also, it is well accepted that the massive galaxies formed their stars at earlier epochs as they appear 
redder at present times, compared to the younger galaxies in which star formation is still going on (Fontanot ct al. 
2009, and references therein). This phenomenon, which are commonly referred as 'downsizing', indicate a possible 
role the massive black holes would have played in quenching the star formation in massive galaxies at earlier epochs 
(Somerville et al. 2004; Granato et al. 2004; Croton et al. 2006). We follow the semi-analytical scheme proposed in 
Granato et al. (2004) for which an approximate analysis is provided in the appendix A of Shankar et al. (2006), using 
which, we can write the rate of change of cool gas in halos as, 

M,oid{t) = -M.nfaiiit) - M.{t) +nfM4t) - CAU{t) , (33) 

where TZf is the return fraction of stars whose value is 0.3 for a Salpetcr IMF. is the mass loss rate, Minfaii{t) 

is the mass infall rate. This yields a stellar mass at present epoch as (see eqn 17 in Shankar et al. 2006), 

= (34) 

where fcosm = 1/6 is the cosmic baryon ratio, fsurv is the fraction of stars surviving up to now and its value is 
approximately 0.6 for a Salpeter IMF. As fsurv, fcosm and TZf are constants, therefore we get Mstar oc Mh/ C. In 

case of SNe feedback the feedback factor is written as £ = °'^^'^e^ ^ ' where Ehind = [v1f{c){l + fcosm) with 
/(c) as given in Mo ct al. (1998), represent the binding energy of gas in the halo per unit mass. The quantity in the 
numerator is nothing but Esn / M — 2v1 where is the value of velocity at the critical point for the case of only SNe 
injection and no AGN. When the AGN is also present we can use the modified value of velocity at the critical point 
given in equation 27 which is v^^-^^ = tj^ — (1 — V)v1/2. Hence in the case of both SNe and AGN feedback we can write 
the loss term C as, 



i«2/(c)(l + /,„,„) \vlf{c){l + f,osm) 



(35) 



Wcrit is the velocity at the critical point and Vc the virial velocity of the halo. While the velocity at the critical point 
measures the strength of SNe and/or AGN as a mechanism of mass expulsion, the velocity at the virial radius is a 
measure of binding energy of the halo. 

For NFW halo, Vc oc M^/^ . Also v1 cx A/, cx A/^-^^ (see equation 32). Using these we get 

°^ ^ 17 mm (36 

l-Tlf+ ./"a Mbreak + M^-''' 

where C, D and Mbreak are constants. From the above relation we can see that for small halo masses we have 
cx M^J'^ and for large halo masses, when the AGN dominates, we get oc M^-^^. The break occurs and the 
SHMR peaks roughly in the range ^ 10^^~^^'^Mq as the ratio I'crit/'^c becomes minimum around this mass. 

In Figure 4, we plot the SHMR at present day, against the halo mass obtained by using equation 34, 35 and 32. The 
set of upper three blue lines correspond to quiescent mode (w^ = 180 km s~^) and the lower three red line represent 
the starburst mode (w* = 500 km s~^). In each set the dashed, solid and the dash-dotted lines correspond to three 
different redshifts of collapse, Zvir = 2.4 and 6 respectively. We stress here that the lines denote the value of SHMR 
at present day (2 = 0) for the galaxies which got viriallized at a particular Zvir. 

The squares, dots, and triangles in Figure 4 represent the stellar and halo mass data inferred from works on TuUy 
fisher relation by Gcha ct al. (2006), Springob et al. (2005) and Pizagno et al. (2007) respectively. The stars and 
diamonds represent the data from (More et al. 2011) and (Conroy et al. 2007) who deduced the stellar and halo mass 
using stellar dynamics. Red circles shows the estimates based on weak lensing by (Mandelbaum et al. 2006). For 
details on the data sets the reader is referred to Leauthaud et al. (2012) and Blanton et al. (2008) and the original 
papers for the the data sets. We have also shown the SHMR obtained by the technique of halo abundance matching 
(Moster et al. 2010; Behroozi et al. 2010; Leauthaud et al. 2012) using broken thick gray lines. At the higher mass 
end we note that our results agree with these works. To be more specific we find a stellar to halo mass slope of 0.26 
as mentioned above, which is in agreement with the value 0.29 deduced in Behroozi et al. (2010). The slope at higher 
mass end depends on the M, — Mh relation, which is still being debated in the literature. However, even if we use a 
different scaling like M, oc cr'' °^ given by Tremaine et al. (2002), then the slope at high mass end becomes 0.33 which 
is also in agreement with observations and other works. 

The slope at low mass end as found by halo abundance matching is roughly ^ 2.2, which is larger than 5/3 from 
theoretical considerations in this work and others (e.g. Dekel & Woo 2003). however, it is possible that other physical 
processes not considered in our model can explain the discrepancy. For example, if v^, depends on halo mass, with the 
efficiency of SNe energy injection being larger in low mass halos, the slopes can be reconciled. 
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Fig. 4. — SHMR at present time from theoretical considerations in this work is compared with the observational data for individual 
galaxies courtesy Alexic Leauthaud. The upper three (blue) lines represent the quiescent mode (d, = 180 km s~^) and the lower three 
(red) lines represent the starburst mode (t>* = 500 km s~^) of galaxies. In these the dashed, solid and dash-dotted lines correspond to 
collapse redshifts, Zvir = 2, 4 and 6 respectively. The data points are from Mandelbaum et al. 2006 : filled red circles, Geha et al. 2006:filled 
green squares, Pizagno et al. 2007 : filled purple triangles, Springob et al. 2005 : dots. More et al. 2011 ; stars, Conroy et al. 2007 : 
diamonds. The thick gray dashed, dotted and dash-dotted lines represent the results of halo abundance matching from Leauthaud et al. 
(2012), Behroozi et al. (2010) and Moster et al. (2010) respectively. 

If one naively compares the lines from our analytical calculation with the observational data points for individual 
galaxies then the plot seems to say something interesting. By looking at the data for low to intermediate mass galaxies 
(Geha et al. 2006; Springob et al. 2005; Pizagno et al. 2007) one may infer that the data points lie systematically below 
the results of halo abundance matching. There are a lot of galaxies which have lower stellar content than expected from 
halo abundance matching. However, if we compare these with our lower three red lines which are for starburst mode 
then these data points can be reconciled. As we have already described, the outflow activity in low to intermediate 
mass galaxies is governed by SNe and starbursts. The two modes we are following in this paper represent two extreme 
efficiencies of energy injection by SNe. Quiescent galaxies like ours lie on one extreme and the violent galaxy like 
M82 on the other extreme. All the galaxies in low to intermediate halo mass range (< 10^^'^Mq) are covered between 
these two extremes. Interestingly in Figure 4 also, most of the data points lie between the blue and red lines of our 
theoretical model which represent the quiescent and starburst mode of energy injection respectively. This is not mere 
a coincidence given the simplicity of our model and shows the importance of outflows in shaping the galaxies. 

In the above analysis, we have assumed that the relation given in equation 34, which is valid in the long time limit, 
gives the stellar mass in the present day universe. In the hierarchical structure formation scenario where small halos 
form first and undergo periods of high merger rates, this assumption may not be completely valid. However, recent 
observations and theoretical works have shown that cosmic downsizing mitigates the effects of hierarchical structure 
formation models, and that massive galaxies are believed to form stars at high redshift after which they evolve passively. 
In this regard, the consistency of our results with observations of massive galaxies becomes important. 

6. DISCUSSIONS 

In this paper we have studied galactic outflows driven by SNe injection and AGN radiation. We have calculated the 
outflow properties in halos ranging from low to high mass. The treatment is analytical which enabled us to extract 
fundamental results for these feedback processes. Here we recapitulate the results and the caveats. 

We find that the NFW gravity causes the outflows to stop inside virial radius in intermediate to high mass halos. 
In Appendix C, we have shown the solutions for winds from NFW halos in terms of Mach number, where the closed 
contours show the importance of NFW halo. In other words, NFW gravity causes the flow to stop inside the virial 
radius of the halo if the energy injection from SNe is not large. We note that a minimum value of SNe energy injection 
so as to produce Vi, Vs is required for the gas to escape the galaxy. Our models in which the gas can not escape the 
galaxy, provides a natural explanation for the circumgalactic material observed inside the halos and also predict mild 
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winds in quiescent star forming galaxies such as our Galaxy. 

When an AGN component is included, the contours in the Mach number-radius plot can open up for massive galaxies 
(Mji > IO^^'^Mq). This implies that AGN radiation can become important in winds from massive galaxies (such as 
ULIGs), reaching a wind velocity of i^wind ^ 1000 km s~^. If we consider that neutral clouds arc entrained in this wind 
then the speeds of the cold clouds can be at the most equal to Wwind- This result is consistent with observed outflow 
speeds in post-starburst galaxies, ULIRGs and Low-ionization quasars at intermediate to high rcdshifts (e.g. Tremonti 
et al. 2007). We have derived a general expression for the terminal speed of the wind from the NFW halo which can 
be written as Wqo = (4wJ + 3(r — l)vl — 4w^)^/^, where is the contribution from SNe and the term with w, stands for 
the momentum injection from AGN. According to our results 1000 km s~"'^is an upper limit for the starburst driven 
winds while in case of AGN driving which is possible in high mass galaxies only, the velocities always exceed 1000 km 
s-i. 

We have also shown that our results can explain the observed trends of stellar to halo mass ratio. We find that the 
stellar mass scales with the halo mass as Af* cx M°-^^ at the higher mass end and Mi, (x M^-^^^ at the lower mass end. 
The slope at higher mass end agrees with the observations and abundance matching results. We find that the large 
scatter in observational data at the lower mass end is due to the fact that the efficiency of energy injection is different 
in different galaxies. We would like to mention here that we have used a simple recipe to calculate the stellar mass 
corresponding to a particular halo mass using the outputs from our wind models. Implementing feedback recipes from 
our wind models into a full semi analytical work is beyond the scope of this paper. 

We note that the reprocessed IR radiation is sufficient to drive strong outflows from high mass galaxies and the 
resulting feedback is enough to explain the mismatch between stellar mass and halo mass at high mass end. This is 
important as one need not to rely on UV radiation which is supposed to be attenuated quickly within a small distance 
from the center. 

The injection of energy and mass in our model occurs in the central region. However recent observations show the 
evidence of outflows emerging from individual star clusters which may be situated away from the center of the galaxies 
(e.g. Schwartz et al. 2006). It becomes however analytically complicated to take into account the contribution from 
these clusters and combine it with the effects of a central AGN. It is however clear that effects of the feedback would 
be higher if mass and energy injection occurs even beyond the central region and from a large number of clusters. 

The coupling between gas and dust is assumed to be strong in our work which is reasonable considering the ionized 
nature of these hot winds. However one may question the survival of dust grains. For this we refer to the observational 
evidence for dust in the spectra of AGNs, which shows that the dust does survive and it may do so by residing in small 
clumps around AGN as discussed in Krolik & Begclman (1988). 

We have considered momentum transfer from AGN radiation in optically thin limit. In actual practice the situation 
can be quite complicated and it may not be exactly accurate to use a constant F. We hope that the simulations of 
winds with full radiation transfer covering from small to large length scales will be able to verify the simple ideas 
presented here. We have used a relation between the black hole mass and the dark matter halo mass. We have 
projected this relation backward in time using scalings of NFW parameters with redshift. There are observational 
evidence that black hole masses at high redshifts are generally higher, however it is hard to predict the correct scaling 
of black hole mass with redshift as the formation and growth of black holes is a complicated problem in itself. 

We have not considered the effects of radiative cooling in the present work. Radiative cooling is important for the 
thermodynamics of the outflow. The GG85 solution with implementation of cooling has been studied in the past (Silich 
et al. 2004; Tenorio-Tagle et al. 2007; Wiinsch et al. 2011; Silich et al. 2011). These calculations conclude that cooling 
(if not catastrophic) does not affect the velocity and density but causes the temperature to decay more rapidly. Wind 
solutions with cooling can not be worked out analytically for general values of parameters and one has to use the 
numerical computation, for particular cases. We have discussed the hot phase of the outflows which is tractable using 
analytic hydrodynamics. However in actual practice the physics of hot, cold and molecular phase of galactic outflows 
might be tangled to each other . The study of all these components demands state of the art numerical simulations. 
We hope to study the multi-phase character of these winds in a future paper. 

To summarize, we have found analytical solutions for SNe and AGN driven winds from realistic dark matter halo. 
Our results show that the two feedback processes operate effectively at two ends of the galaxy luminosity function as 
expected. The wind velocities for escaping winds resulting from our calculations explains a variety of observations. We 
flnd that AGN can drive the gas to speeds > 1000 km s~^. We find an intermediate mass range in which the outflows 
can be highly suppressed and for these halo masses the gas can not escape into the IGM. However, the gas is still 
driven to large distances within the halo. This result explains the resent observation of gas reservoirs in our Galaxy 
and other galaxies. Using the results of our analytical models, we have derived a stellar mass to halo mass relation 
using simple recipes. We find the derived M^^-Mh relation matches the results of observations and halo abundance 
matching. Thus, our findings provide a possible explanation for missing baryons in galaxies. 

We thank Mitchell Begclman for critical reading of the manuscript and Rachel Somerville for useful discussions. We 
are grateful to Alexie Lcauthaud for providing the observational data used in Figure 4. 
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APPENDIX 

APPENDIX A: DETAILED DERIVATION OF WIND EQUATION 
Substitution of dp/ dr from the continuity equation 1 into the momentum equation 2 yields 
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from energy equation we get 



7 dr 

dv 
jv dr 

dv 



7 XT 

1 del 
7 dr 

1 del 



1 dv 
V dr 



dr ^ — 1 dr 



7 
pv 



m \ 


mv 


-fir) 


d^{r) 


pv) 


P 


dr 


rh ^ 


mv 




d$(r) 


pvj 


P 


dr 



After a httle manipulation we obtain 
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From the above two equations we get the following wind equation, 
M"^ - 1 dM^ 
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where e(r) = ^ + 7^^- This is the general form of wind equation with contributions from energy and mass injection, 

gravity and external driving force. For the solar wind case, $ = — GMQ/r and the corresponding analysis can be found 
in Lamers & Cassinelli (1999). 

APPENDIX B: SUBSONIC PART OF SNe & AGN DRIVEN WIND 

Here we show the subsonic part of the solution for SNe and AGN driven wind. Apart from the energy and mass 
injection we take /(r) = FGM./r and <J>(r) = $,(r) = —GM,/r. Substituting these values in equation 4 or A8 we 
obtain the following wind equation. 

From continuity equation we obtain pvr"^ = mr^/3. If we substitute this in energy equation we get, 
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In the subsonic regime r « r^. Therefore in this case the third term on left hand side which is due to NFW gravity 
can be neglected as shown below, 
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Neglecting the third term, (^ihil d^MiM. 



ln(l + r/r,) 



3 

■mB? 



3 [ln(l + c) - c/(l + c)] V 1 + rjrs 
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where we have defined v, = \J GM,/2R, a velocity characteristic of the black hole. In the above equation for r << i? 
we can neglect all the terms with or and then dividing by rhR^ /3 we obtain. 
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substituting this value of e(-R) in equation B3, we get 
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Substituting this e(r) and also /(?') 
can get. 



TGM,/r^ , $ = $, + $NFW, 7 = 5/3 in wind equation 4 of the main text, one 
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Solution of this equation gives the mach number and velocity in the subsonic part of the wind. We can clearly see from 
the above equation that for v, ~ which will mean the AGN is not present, the subsonic part of the CC85 solution 
is recovered. Also the same thing happens when F = 1 because in that case, although the AGN is present but the 
outward radiation force cancels the inward gravitational force everywhere. 

APPENDIX C: MACH NUMBER VERSUS DISTANCE DIAGRAMS 
Integral of the wind equation 26 for the winds with NFW gravity and AGN is, 
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2v1 — (1 — r)i'^. The above equation can equivalently be written as 
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where we have used the definition of function 5y[M) given in §2.1. In this equation A is an arbitrary constant. If 
we set A = l/(2i?iJcrit)^, then for r = i? and M. ~ \, both LHS and RHS of the above equation become equal to 1. 
In Figure 5 we plot the contours of Mach number versus the radius for three galaxies with a different halo masses. 
The upper three panels represent the winds without AGN (w, = 0) and the lower three show the effect of the AGN 
momentum injection. Different contours in each panel correspond to a different value of the constant A. The thick 
blue contour in each panel represent the case with a critical point at r = i? = 200 pc and this one is used in the main 
text to calculate wind properties. 

With the inclusion of NFW gravity one encounters a wall type of behaviour at a particular r. Beyond the wall the 
solutions real and physical are not possible. We would like to mention here that similar situation arises in adiabatic 
solar wind problem as shown in Panel (c). Figure 2 of Holzer & Axford (1970). The difference is that in the case of 
galaxy, the energy injection causes a critical point at R=200 pc. From Figure 5 one can infer the interesting fact that 
the AGN is not able to drive the gas out of galaxy for intermediate halo masses but it can do so in high mass galaxies. 
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Fig. 5. — Mach number versus the radial distance for the winds from a NFW dark matter halo. Upper row of plots shows the effect of 
NFW gravity only and the lower row takes into account both the NFW gravity and an AGN. The three plots in each row are for three 
different halo masses, x-axis in each plot extends upto virial radius corresponding to the halo mass used. The thick blue line represents 
the solution with a critical point at R = 200 pc, which is used in this work. We have taken ii* = 180 km . 



APPENDIX D: POSSIBLE EFFECTS OF RADIATIVE COOLING 

In this section we discuss the effect of radiation loss on the wind properties. So far we have neglected the effect 
of radiative cooling, since it is known that the energy loss due to radiation is too small to affect the dynamics of 
winds (Grimes et al. 2009). However the cooling does effect the thermodynamics of the wind and can cause steeper 
fall in temperature as compared to the case with pure adiabatic expansion. We show this here with a simple case of 
supersonic escaping wind without gravity for which the wind equation with radiative cooling can be written as, 



7W2(7V{2(2/3) + 2) dr r 2pv e{r) 

^ r ^ IGtt e(r)2A12 



(Dl) 



where we have used p = M/{ATTvr'^), 7 = 5/3 and — 2e{r)/{Ai'^ + 3). Here A(r) is the cooling rate in erg cm'^ 
s~^. We use A oc T"-^ for the temperature range 10*^ < T < 10^ and A cx in the temperature range 10"* < T < 10^ 
(Sutherland & Dopita 1993). To solve this equation we have to supply e(r) and M. The specific enthalpy, e(r) can be 
obtained from the integral of energy equation which is given below, 

eir) = '- + ^ = ^-1 r ^(W)dr = ^ - (D2) 

Considering the fact that total energy extracted by cooling over the entire path is ~ 10% of the adiabatic losses (Wang 
1995; Grimes et al. 2009) we can neglect Ecooi as compared to Esn- Thus we get e(r) = Esn/M — 2v1. Using the 
e(r) and a value for M in equation Dl, we can solve it to obtain mach number as a function of distance. 

In Figure 6 we plot the sound speed Cg = [2e(r)/(A^2 _|_ 3)]^/^ against r, where M. is obtained by solving equation 
Dl. We have solved equation Dl for the quiescent mode (w^ = 180 km s^^) and the starburst mode (w* = 500 km 
s"^). We have used M ~ 1 MQ/yr and 10 Mo/yr as two fiducial value of mass loss rate for the quiescent mode and 
the starburst mode respectively. The sound speed corresponding to quiescent mode is shown by dotted line and that 
of the starburst mode by a dashed line. We have also shown the corresponding cases without cooling using a solid 
line and a dash-dotted line. By comparing the sound speeds with and without cooling, we find that the cooling causes 
the temperature to decay more rapidly. Also cooling is effective at smaller distances (r < 10 kpc) as the density and 
temperature are larger there. Wind speed is given by the relation e(r) = v'^/2 + Cg/(7 — 1) ~ 2v'^. The terminal 
wind speed is obtained by neglecting the sound speed which is known to decrease with distance. As the sound speed 
go down even more rapidly in case of cooling, therefore the wind speed beyond a distance of 10 kpc does not differ 
from the case without cooling and is given by, Wwind ~ 2Vi,. We would like to mention here that qualitatively similar 
and quantitatively more accurate result has also been obtained by numerically solving the basic fluid equations with 
cooling for winds from individual star clusters in Silich et al. (2004) and Tenorio-Tagle et al. (2007). Whether the flow 
undergoing radiative energy losses can achieve a steady state, is also an interesting problem. We refer the reader to 
Silich et al. (2003) where the time dependent problem on a 2-D grid has been attempted and the result show that the 
flow can become steady after sufficient amount of time. 
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Fig. 6. — Sound speed for the winds with cooling arc compared with the cases without cooHng. The Solid and dash dotted lines represent 
the case without cooling for quiescent mode and starburst mode respectively. Their counterparts with cooling are shown by dashed and 
dotted lines. 
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